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Cold dust grains are responsible for the far infrared and submillimetre (FIR/submm) emission observed by Herschel and Planck. 
CO ' Their thermal emission is usually expressed as a modified black body law in which the FIR/submm dust opacity, or mass absorption 

coefficient (MAC), is described by the mass absorption coefficient at a given wavelength, K^g, and the temperature- and wavelength- 
independent emissivity spectral index, p. However, numerous data from previous space and balloon-born missions and recently from 
' J ' , Herschel and Planck show that the dust emission is not well understood, as revealed for example by the observed anti-correlation of 

■ j8 with the grain temperature. 

r ^ ' The aim of this work is to measure the optical properties of interstellar dust analogues at low temperatures to give astronomers the 

I necessary data to interpret FIR/submm observations such as those from the Herschel and Planck satellites. 

. We synthesised, via sol-gel methods, analogues of interstellar amorphous silicate grains, rich in Mg and Ca, and having stiochiometry 

of olivine and pyroxene. The samples are characterised by various techniques to determine their composition, size, amorphization 
degree. All the amorphous samples are annealed at 1 100°C to study the crystallised materials for comparison. We measured the MAC 
of all the samples, in the 2-25 yum range at room temperature and in the 100-1000/1500 yum range for grain temperatures varying from 
300 to 10 K. 

^ ■ The experimental results show that for all the amorphous samples, the grain MAC decreases when the grain temperature decreases and 

^ ' that the local spectral index, /?, defined as the slope of the MAC curve, is anti-correlated with the grain temperature. These variations, 

which are not observed in the crystallised samples, are related to the amorphous nature of the samples. In addition, the spectral shape 
of the MAC is complex: it cannot be described by a single spectral index over the 100-1500 yum range. At short wavelengths (A < 
^ ■ 500/700 /jm),/? is in the range 1.6-2.1 for all grain temperature and grain composition. However, at longer wavelengths (A > 500/700 

QQ ' yUm), /? < 2 for samples with a pyroxene stoichiometry and yS > 2 for samples with an olivine stoichiometry. 

1^-^ I The dust properties in the FIR/submm domain and at low temperature are more complicated than expected. The simplifying asymptotic 

, expression based on a single temperature- and wavelength-independent spectral index used by astronomers is not appropriate to 

■ describe the dust MAC and hence the dust emission, and may induce significant errors on the derived parameters such as the dust mass 
• ' and the dust physical and chemical properties. Instead, dust emission models should use the dust MAC as a function of wavelength 

and temperature. 

o : 

Key words. Astrochemistry - Methods: laboratory - Techniques: spectroscopic - (ISM:) dust, extinction - submillimetre: ISM- 
Infrared: ISM 

1. Introduction the silicates in the ISM are completely amorphous dKemper et al.l 

'2004). It is more difficult to study amorphous silicates than 
C3 The far infrared and submillimetre emission of the Interstellar crystalline sihcates because the vibrational bands at ~ 9.7 /jm 
Medium (ISM) is dominated by grains of size of the order of a and ~ 18 yum of amorphous silicates are broad and sti'ucture- 
few tens to a few hundreds of nanometers, mostly composed of jgss whereas in crystalline silicates, they exhibit numerous sharp 
silicates and oxides ([Desert et al. 1990; Compiegne et al. 201 1). spectral features due to the long-range order of the crystalline 
These grains ai-e ubiquitous. Their physico-chemical proper- lattice. Consequendy, the composition of the amorphous sil- 
ties and their evolution within the Ufecycle of the ISM are jcates is not as well constrained as for crystalline silicates, 
constrained by the spectroscopic observations in the mid in- Nevertheless, the comparison of astrophysical specti-a with labo- 
frared spectral range (2-50 ^m) via the study of the vibration ratory spectra of amorphous silicate dust analogues and the con- 
of the chemical bonds constituting the grains. Results from straints coming from the depletion of the elements from the gas- 
the space-born satellites ISO (Infrared Space Observatory) and phase, indicate that the amorphous silicates should have an av- 
Spitzer show that although crystalUne silicates such as forsterite grage composition close to the one of the crystalline silicates: 
(Mg2Si04) and enstatite (MgSiOs) are observed around stars, Mg2Si04, MgSiOj. In addition dust grains certainly contain 
iron, which is highly depleted, but the form in which it is in- 
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in oxides, in the form of iron inclusions). They may also contain 
other minor elements highly depleted from the gas-phase such 
as calcium or aluminium. These elements may be incorporated 
within the silicates chemical network or in the form of oxides. 
Observations in the far infrared and submillimetre range trace 
cold dust (T~10-40 K). It allows us to detect pre-stellar cores, 
prior to their collapse and thus to study the very first steps of star 
formation. It is used to estimate the dust mass and density of in- 
terstellar cold clouds. When the interstellar grains are big enough 
(above 10 nm), the absorption of UV-photons from the ambient 
radiation field does not induce a significant grain temperature 
variation and these cold grains radiate at thermal equilibrium. In 
cold envkonments of the ISM, this emission dominates above 
other dust component emission at wavelength greater than 100 
fj.m. Assuming an optically thin medium, which is usually the 
case in the FlR/mm wavelength range, their modified blackbody 
emission is expressed as: 

1(A) = B(^, T) t(A) = B(^, T) k(A) Md d"^ Q"' (1) 

where 1) is the specific intensity and B )(T) the Planck func- 
tion at the dust temperature in W.m^^.Hz^'.sr"', the mass 
of dust along the line of sight, d the distance of the source, Q 
the beam solid angle and k(A) the mass absorption coefficient 
(MAC, cm^.g"'). Simple semi-classical physical models of ab- 
sorption, provide, at long wavelength, a temperature independent 
asymptotic behaviour for the dust emissivity that is expressed as 
follow: 

K = (7 (2) 

Ao 

where the emissivity spectral index, /3, is ranging from 1 to 2, a 
value of 2 being commonly adopted in astrophysical studies. 
Before the launch of the Herschel and Planck satellites the, rel- 
atively few, observations in the far infrared and submillimetre 
wavelength revealed that the dust emission was not easy to in- 
terpret. Among the observed anomalous behaviours are the far- 
infrared excess observed by the instrument FIRAS on board 
the COBE (Cosmic Background Explorer) satellite (Reach et al. 
1995) and the /3-T anti-correlation derived from the observations 
of the two balloon-born missions PRO NAOS and ARCHEOPS 
(iDupac et al.ll2003l:lDesert et al.ll2008h . The FIRAS observations 
revealed the existence of a significant millimetre excess in the 
dust emission spectrum, with respect to a modified black body 
law with a spectral index /3 equals to 2. The data can be well- 
fitted with a spectral index value of 1 instead of 2. It was pro- 
posed that this excess is due to the presence of very cold dust 
(~ 9 K), in addition to a warmer compo nent responsible fo r 
the emission maximum near 100 fim (Finkbe iner et al.l ri999h . 
However, the fact that the excess is well correlated with the FIR 
dust emission indicates that it could originate from unidentified 
processes intrinsic to the grains radiating in thermal equilibrium, 
without requiring an additional colder component. The deviation 
of the spectral index value from the classical value of 2 was also 
observed in various interstellar regions for which the p value 
increases (up to 3 or 4) when the mean dust temperature de- 
creases (down to 10 K) whereas it is as low as 1 when the dust 
temperature equals 80 K (Dupac et al. 2003; Desert et al. 2008). 
The anti-correlation of /3 and T is now confirmed by the first re- 
sults of the Herschel mission which show that the value of the 
spectral index, derived from the fit of the observed spectral en- 
ergy distribution (SED) of regions in our Galaxy and in exter- 
nal galaxies, often deviates from the value of 2 (Anderson et al. 
l20Ta iRodon et alJ 12011 lBot etal.1 lloioh . The Planck Early 
Results confirm the dust emissivity variations and indicate that 



they are ubiquitous. First results show that the dust emissiv- 
ity index in the solar neighbourhood of our Galaxy changes 
from y6-1.8 in the FIR to a flatter beh aviour in the submillime- 
tre ( Planck Collaboration et alJl201 Id) . They also clearly show 
that dust in nearby galaxies (the Large and Small Magellanic 
Clouds, LMC and SMC in particular) has a significantly flat- 
ter FIR emissivity {P-1.5 and 1.2 respectively) than our Galaxy 
(Hangk Collaborati on et al. 201 lb). Complementarily, using the 
new Herschel data. lParadis et alj (i2010) showed that the emissiv- 
ity index of dust in our Galaxy actually varies with temperature 
in the far infrared. With its full sky survey and unprecedented 
sensitivity, Planck has also allowed us to detect more than 10 000 
cold clumps in the Galaxy (Planck Collaboration et aL 2011. J)- 
The temperature of the Planck cold clumps spans from 7 K to 
17 K, with /3 ranging from 2.8 to 1.4. It has been shown that the 
data are not consistent with a constant value of the dust spec- 
tral index /3 over the whole range of temperature, confirming 
the anti-correlated behaviour observed previously. Furthermore, 
Monte-Carlo simulations have demonstrated that such an anti- 
correlation cannot be explained by a fitting degeneracy given 
the error bars, and could be a real feature of dust emission 
dPesert et alJl2008t IVeneziani et ani2010t iParadis et alJl2010l) . 
Motivated by the results from infrared satellites, many exper- 
imental and theoretical studies have been undertaken to char- 
acterise the optical properties of silicate dust analogues in 
this spectral range. Most of t h e exp erimental studies focus 
on the MIR range dJager et alj d 1994b and associated papers, 
iKoike et al.l (l200Cl) ) alfliough in some cases, experimental spec- 
tra were reco rded up to 100 or 300 //m and sometimes at low 
tem perature s (Henning & Mutschke' 1997: Chi hara et alj|200Tl: 
Koike et al.ll200 6: Mutschke et al. 2008). However only a few 
studies were dedicated to the far infrared and submillimetre 
spectral range. In the far infrared range, the mass absorption 
coefficient (MAC) may be influenced by various factors such 
as the dust aggregatio n, the grain size , the p r esence of con- 
ductin g particles (see iStognienko et al.l (Il995h : iHenning et al.l 
(Il995h for a modelisation of these effects). The dust temper- 
ature may also play an important role as demonstrated by the 
experimental studies from Agladze et al. ( 1996); Mennell a et al.l 
( 1998); Boud et et al.l (120051) . I Agladze et alJ (11996) measured the 
absorption spectra of crystalline enstatite and forsterite and of 
amorphous Mg-rich submicronic silicate grains at very long 
wavelength (700 /^m -2.9 mm) and very low temperature (1.2- 
30 K). Mennella et al. ( 1998) performed similar measurements 
in the 20-2000 yum range on forsterite (Mg2Si04) and amorphous 
and cr ystalline fayahte (Fe 7Si04) at higher temperatures (24- 
300 K). lBoudet et alj(l2005h studied amorphous MgSiOs and sil- 
ica (SiOa) submicronic grains in the 100-2000 yum range and low 
temperature (10-300 K). These three studies all show that the 
absorption coefficients of the amorphous samples are frequency- 
dependent and temperature-dependent. However, because these 
studies were made on different materials, synthesised with dif- 
ferent techniques and characterised at different temperatures and 
spectral ranges, they are difficult to compare. 
We present in this paper an experimental study of the optical 
properties of interstellar silicate dust analogues in the wave- 
length range 100-1500 fim and from 300 to 10 K. The studied 
samples have a mean composition close to Mg2Si04, MgSiOs 
and CaMgSi206. We study amorphous and crystallised samples 
to investigate the influence of disorder on the variation of the op- 
tical properties of the material as a function of the grain temper- 
ature. The synthesis of the analogues and their characterisation 
are presented in Sect. |2] The experimental details related to the 
spectroscopic measurements are described in Sect. [3] The results 
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Fig. 1. TEM bright field images of (a) sample Fl, (b) sample F2 and (c) sample E. All samples are prepared by a sol gel method. 
The microstructure consist of small clusters stick together. The sample have a high porosity and a high surface/volume ratio. 

Table 1. Characteristics of the synthesised dust analogues. 



sample 


sol-gel synthesis 


composition 


specific area 


(p pores 




recipe 




K.g) 


(nm) 


Fl 


TEOS + Mg(0Et)2 + H2O + HCl 


Mg2 3Si04 




-5-10'" 


F2 


TEOS + Mg(0Et)2 + H2O + C4H5NH2 


Mg28Si04 




~5-10(" 


F3 


TEOS + Mg(N03)2 + H2O + HNO3 


Mg2.05SiO4 




~ 5-100 


E 


TEOS + Mg(N03)2 + H2O + HNO3 


Mgo.wSiOs 


132 


9(2) 


D 


TEOS + Mg(N03)2 + Ca(N03)2 + H2O + HNO3 


Cao.98Mgo.9Si206 


83 


17(2) 



Notes. ''^ Estimated from TEM images. *'* Derived from BET measurements. 



are presented in Sect. |4] and discussed in Sect.|5] The astrophys- 
ical implications are discussed in Sect. |6] 



2. Sample synthesis and characterisation 

2.1. Sample synthesis 

We have chosen to study different amorphous Mg-rich silicates 
representative of astronomical silicate dust (Table [T]i: 3 samples 
with a stoichoimetry close to olivine (Mg2Si04) and 2 samples 
with a stoichiometry close to pyroxene ((Mg,Ca)2Si206). The 3 
silicates with an olivine stoichiometry were synthesised by sol- 
gel method using three different procedures. The sol-gel process 
is based on the hydrolysis and condensation of tetraethoxysi- 
lane (TEOS) in a solvent using water and a Mg-rich compound 
as the reactant, in presence of a catalyst. The main differences 
between the 3 adopted procedures are the choice of the Mg- 
rich compound, which can be a nitrate salt (Mg(N03)2) or an 
aUcoxide (Mg(OEt)2), and the nature of the catalyst which can 
be acidic or basic. The first sam ple, called Fl , was synthe- 
sised following the procedure of iMitchell et al.l (Il998h using 
Mg(OEt)2 and HCl. The second sample, called F2, was synthe- 
sised following a procedure from Coupeaud (unpublished) using 
Mg(OEt)2 and butylamine (CH3(CH2)3NH2, hereafter BA). For 
this method, a solution of TEOS (3 mmol, 0.12 M) in tetrahy- 
drofurane and 2-methoxyethanol was prepared at room tem- 
perature with a TEOS:BA:H20 ratio of 1:0.5:4 (lEl Hawi et al.1 
[2009). A solution of Mg(OEt)2 (6 mmol, 0.24 M) prepared as 
described bv IMitchell et all (1 19981) was then added to the TEOS 
solution and the mixture was stirred and heated at 50°C for 
12 days. Finally, the obtained suspension was centrifuged and 
washed with tetrahydrofurane and 2-methoxyethanol. The third 
sample, cal led F3, was synthesised following the procedure of 
iGillot et all (|2009) using Mg(N03)2 and HNO3. The samples 
with a pyroxene composition, called E and D, were synthesised 



using the sol-gel method from iGillot et al.1 (l2009l) using TEOS 
and Mg(N03)2 with a ratio of 1 for the E sample and using 
TEOS, Mg(N03)2 and Ca(N03)2 for the D sample. To remove 
the organics from the synthesis, all the samples were heated 
up to 400°C or 500°C depending on the synthesis procedure. 
The composition of the synthesised samples were determined 
by Inductively Coupled Plasma Mass Spectroscopy (ICP-MS, 
Antellis Company), they are indicated Table [T] All the samples 
contains OH groups and water coming from the sol-gel method 
that partly remain within the samples after thermal annealing 
at 400-500°C. The influence of OH groups and water on the 
FIR optical properties of the samples is discussed in Sect. |5] In 
addition, about half of each of theses amorphous samples was 
annealed for 24 hours at 1100°C in oven under ambient atmo- 
sphere. The resulting samples, called Flc, F2c, F3c, Ec, and Dc, 
are crystalline as shown by the XRD and NMR analysis and MIR 
spectroscopy (see Sect. 12.41 Sect. l2.5l and Sect. 14. lb . Thus, a total 
of 10 samples are considered in this study: the 5 amorphous sam- 
ples produced by sol-gel method and annealed at 400-500°C and 
the 5 crystallised samples obtained by annealing of the amor- 
phous samples at 1100°C. 



2.2. Specific area and porosity 

The surface area and average pore size of the E and D samples 
were determined by physical gas adsorption technique. An inert 
gas, nitrogen, is adsorbed on the surface of a solid material (on 
the outer surface and, in case of porous materials, also on the sur- 
face of the pores) at a temperature of 77 K leading to a so-called 
adsorption isotherm. It is standard to apply Brunauer-Emmett- 
Teller (BET) method to derive the surface area and pore size 
from the adsorption isotherm. The results show that the synthe- 
sised sol-gel silicates E and D are porous, with an average pores 
diameter between 10-20 nm, and present a high specific area in 
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Table 2. Composition of the crystalline phases measured by 
XRD after annealing at 1 100°C of the silicate samples. 



phase 


phase 


Flc 


F2c 


F3c 


Ec 


Dc 


name 


composition 


phase abundance (in mass %) 


periclase 


MgO 


38 


18 


6 






forsterite 


Mg2Si04 


38 


82 


94 


15 




orthoenstatite 


MgSiOj 


24 






64 




clinoenstatite 


MgSiOa 








21 




merwinite 


Ca^MgSizOg 










7 


diopside 


CaMgSijOfi 










65 


akermanite 


Ca^MgSizO, 










28 



the range 80 - 130 m^/g. For the Fl, F2 and F3 samples, no BET 
measurements were performed. 

2.3. Transmission electron microscopy 

Transmission electron microscopy (TEM) is used to investigate 
the samples prepared by slow evaporation of a drop of the col- 
loidal solution deposited onto a carbon covered copper grid. 
The TEM analysis were performed at the "Service Commun 
de Microscopic Electronique de I'Universite Paul Sabatier" 
(TEMSCAN) on a JEOL JEM 1011 electron microscope op- 
erating at 100 kV with a point resolution of 4.5 A and at the 
"Centre Commun de Microscopic de Lille" on a FEI Tecnai G2- 
20 twin microscope operating at 200 kV. The samples E and D 
(having a pyroxene stoichiometry) have an homogeneous com- 
position. The TEM images of these two samples show that the 
particles present a complex texture with a fine apparent porosity 
with pore sizes of 5-20 nm, compatible with the BET analysis 
(see Sect. 12.2b . The 3 samples having an olivine composition are 
different from each other The Fl sample is chemically inhomo- 
geneous: Mg-rich nanocrystal of ~ 5 nm (presumably MgO) are 
identified on the TEM images together with regions deficient in 
magnesium (Mg/Si ~ 1). The F2 sample is relatively stoichio- 
metric with a slight Mg excess compare to Si. The F3 sample 
is partly crystallised. As for the samples with a pyroxene sto- 
ichiometry, samples with an olivine stoichiometry are charac- 
terised by a fine porosity with an apparent pore size of ~ 5-10 
nm. However no BET measurements were performed on the F 
samples to confirm this size. 



2.4. X-rays diffraction 

X-rays diffraction (XRD) measurements are used to monitor the 
amorphous degree of the samples and to characterise the compo- 
sition of the crystalline phases formed after annealing. The sam- 
ples were ground in an agate mortar and were analysed in the 
Panalytical MFD Pro powder diffractometer at the Laboratoire 
de Chimie de Coordination (CNRS, Toulouse) and in a Seifert 
XRD 3000TT diffractometer at the Centre d' Elaboration de 
Materiaux et d'Etudes Structurales (CNRS, Toulouse). The X- 
ray diffraction measurements were carried out in 0-9 configura- 
tion and diffraction profiles were obtained using a step interval of 
respectively 0.017 ° and 0.015 ° in 26 according to the diffrac- 
tometer used. The analysis show that the samples annealed at 
400 or 500°C are amorphous, except for the F3 sample which 
XRD spectrum exhibits peaks characteristic of forsterite. This 
is confirmed by the NMR measurements and the MIR spectrum 



(see Sections [2 .5 1 and l4m i . 

The quantitative analysis of the phases constituting each crys- 
tallised samples was perfor med by Rietveld refin ement of the 
powder diffraction patterns ( Bish & Howardll988h . The identifi- 
cation of the phases constituting each compound was performed 
using the Eva software (http://www.bmker-axs.de/eva.html}. The 
identified phases were then used as an input for Rietveld refine- 
ment, using the software MAUD (http://www.ing.unitn.it/maud/). 
This technique allows us to derive the relative concentration of 
the different crystalline phases with an error of about 5%. The 
Rietveld refinement of the Flc, F2c, F3c samples XRD mea- 
surements shows the existence of two crystalline phases, peri- 
clase (MgO) and forsterite (Mg2Si04), respectively in the fol- 
lowing ratios : 38%/38%, 18%/82% and 6%/94% (Table |2li. In 
addition, around 24% of orthoenstatite (MgSiOa) is present in 
the Flc sample. The powder pattern of the sample Ec shows that 
the sample is a mixture of forsterite (15%), orthoenstatite (64 
%) and clinoenstatite (21%). The powder pattern of the sample 
Dc shows that the sample is a mixture of 3 crystalline phases: 
merwinite (7%), diopside (65%) and akermanite (28%). 

2.5. Nuclear magnetic resonance 

Solid-state NMR spectra of ^^Si were recorded at the Laboratoire 
de Chimie de Coordination (CNRS, Toulouse) on a Bruker 
Avance 400WB spectrometer equipped with a 9,4 T magnet. 
NMR spectra are referenced to TMS. For all experiments we 
used a 4 mm diameter probe and a magic-angle spinning (MAS) 
rate of 12 kHz. MAS ^^Si NMR spectroscopy using high power 
decoupling (HPDEC) and 'H-^^Si cross polarisation (CP) se- 
quences is a powerful tool to characterise the silicate structure at 
the atomic scale. The chemical shift data of ^^Si NMR spectra of 
silicate samples allow the identification of the different Q"-type 
silicon signals (n represents the number of siloxy groups (-OSi-) 
per silicon atoms, Si-O-Si). Thus, Q"* corresponds to Si(OSi)4, 
Q^^ to -Si(0Si)3, Q- to =Si(0Si)2, Q' to =Si(OSi) and Q" to iso- 
lated tetrahedron Si04. 

The analysis of crystalline samples is compatible with the XRD 
results. The NMR spectra exhibit narrow peaks with chemical 
shifts attributed to Q" and groups that are characteristic of 
respectively forsterite and pyroxene (enstatite and diopside). 
In amorphous samples, these peaks are broader than in crys- 
tals. In addition, because the connectivity in amorphous and 
crystalline materials is different, the amorphous samples present 
a wide distribution of Q" groups. The proportions of the Q" 
species were obtained from the deconvolutio n of the ^^Si 
HPD EC NMR spectrum with the Dmfit program dMassiot et al.l 
|2002). The relative amount of Q" groups are correlated to the 
structure of the sample. Indeed, the F samples are characterised 
by a majority of Q° and Q' moieties (~60%) and exhibit the 
features of Q^, and Q'* too. In addition, the spectrum of F3 
shows a spectral feature of Q*^ group which is very narrow and 
shifted in comparison with the amorphous Q" type silicon signal 
confirming that the sample is partially crystalline and contains 
~14% of forsterite. The deconvolution of the spectrum of the E 
sample reveals a majority of Q' and groups (-- 70%) and the 
presence of and groups. For the D sample, the distribu- 
tion of Q" species is made up of - 40% of Q- and ~ 60% of 
Q^* and Q"^ mixture. In this sample, the lack of Q' group shows 
a high degree of connectivity of the Si04 tetrahedra. These anal- 
ysis thus show that the amorphous samples have a large range 
of connectivity, on average larger than in the associated crys- 
talline structure (Q" for forsterite and for pyroxene enstatite 
and diopside). In addition, the F samples (Mg2-2.8Si04) display 



4 



A. Coupeaud et al.: Low temperature FIR and submm mass absorption coefficient of interstellar silicate dust analogues 



a degree of connectivity lower than the E and D samples because 
the dominant species are Q° and Q' while it is Q' and for the 
E sample, and and for the D sample. 



3. Infrared and submillimetre spectroscopy 

3.1. Sample preparation 

To perform the measurements in the MIR spectral range, we 
mixed our silicate powder with KBr powder (Aldrich) in a mor- 
tar and we applied on this homogeneous mixture a pressure of 
7 tons during 10 min, thus obtaining a 10 mm diameter pellet. 
Typically, we use 0.2-0.5 mg of samples embedded into 100 mg 
of KBr. 

For the FIR spectroscopic measurements, the silicate grains were 
mixed with polyethylene (PE) powder at 120°C and we applied 
on this homogeneous mixture a pressure of 7 tons during 5 min 
to obtained 10 mm diameter pellets. The mass absorption effi- 
ciency (MAC), or opacity, of the dust analogues decreases with 
increasing wavelength. Therefore, in order to cover the spectral 
range from 100 fj.m to more than 1 mm, it is necessary to per- 
form measurements on pellets containing more or less material 
depending on the spectral changes of the MAC. Typically, two or 
three pellets were needed, containing from 10 mg up to 100 mg 
of samples and 100 mg to 250 mg of PE. 

3.2. Experimental set-up 

The MIR/FIR/submm optical properties of the samples were 
measured with the experimental set-up ESPOIRS (Etudes 
Spectroscopiques des Proprietes Optiques InfraRouge et 
Submillimetrique of interstellar dust analogues). It is composed 
of a Bruker IFSl 13 V Fourier Transform InfraRed (FTIR) spec- 
trometer coupled with a continuous flow helium cryostat which 
allows us to cool down the samples from 300 K to 10 K. The 
MIR transmission spectra were recorded using a KBr beamsplit- 
ter, a Globar source and a Deuterated TriGlycine Sulfate (DTGS) 
detector with a KBr window at a spectral resolution of the in- 
strument of 1 cm"'. In the MIR spectral range, the transmission 
measurements were performed only at ambient temperature for 
all the samples. 

To perform measurements in the FIR/submm spectral range, 
the spectrometer is equipped with a Hg lamp source, different 
Mylar beamsplitters with thickness of 23, 50, 125 fj.m and a 4 K 
Ge-bolometer detector. This configuration allows us to perform 
measurements in the 100 - 800 /urn wavelength range (5000 - 
12.5 cm"'). To complete our wavelength range above 800 yum 
we performed experiments on the AILES (Advanced Infrared 
Line Exploited for Spectroscopy) beamline at the synchrotron 
SOLEIL (Brubach et al...201 0). The experimental set-up is com- 
posed of a FTIR spectrometer Bruker IFSl 25 coupled with a a 
Helium pulse tube cryostat (Cryomech PT405). Thanks to the 
use of a 1.6 K bolometer and to the high brightness of the syn- 
chrotron IR light (100-1000 times greater than thermal source) 
we extended the spectral coverage to the lowest energy range, up 
to 1500 yum (6.6 cm"'), with a very good signal/noise ratio. The 
sample holder in the two cryostats has two positions allowing us 
to record the transmission spectra of the sample and the blank PE 
pellet during the same cooling cycle. The transmittance spectra 
were recorded at 300, 200, 100, 30, and 10 K. 



3.3. Data analysis 

The transmittance spectra are calculated by dividing the spec- 
trum of the sample pellet by the one of a blank pellet containing 
the same mass of matrix (KBr or PE depending on the spec- 
tral range) as the sample pellet. Hence, the blank and sample 
pellets do not have the same total mass nor the same thickness. 
The transmission of PE does not vary with the temperature from 
300 K down to 10 K except around -125 fim where PE ex- 
hibits an absorption band which shape is temperature depen- 
dent. The transmittance spectra are affected by two problems. 
Firstly, they present very noisy regions that correspond to null 
absorption of the 23 fim beamsplitter around ~ 155-180 fim. In 
this spectral region, and around ~ 125 fim, we replaced the ex- 
perimental data by a linear interpolation. Secondly, the multiple 
reflexions of the incident FIR beam between the parallel faces 
of the PE-pellets cause interference fringes that are most im- 
portant at long wavelength (typically above 160 yum) where the 
material absorbs weakly. We applied two different methods to 
re move these fringes. T he first method is similar to the one used 
bv lBoudet et al.l (l2005h and consists in performing a polynomial 
fit of the transmittance spectra. We used the following function: 
f{x) - ax^ + bx^ + cx* + dx^ + ex^ + fx + g + h/x. The second 
method is based on the convolution of the transmittance spectra 
with a sine function (C. Tucker, private communication ). 

The MAC of the pellet (sample and matrix), ^-peiiet (cm^.g"'), 
is calculated from the transmittance spectra with the following 
formula (e.g. Bohren & Huffman (1998.) ): 



'^pellet 



M 



X ln(T) 



(3) 



where S is the pellet section (cm^), M the sample mass in the 
pellet (g) and T the transmittance spectrum. We have assumed 
that scattering effects may be neglected since the grains size is 
much smaller than the wavelength. The MAC of the grains, k, is 
related t o the MAC of the pell et by the following relation (Eq. 
(6) from lMenneUa et al.1 (Il998h ): 



k(A, f) = g(e, f) * /fpellets(/l, f) 



(4) 



where f is the volume filling factor of the grains within the matrix 
and g(e, f) is the finite concentration reduction factor The dilu- 
tion factor adopted to make the pellets is high enough to ensure 
that the filling factor f never exce ed 0.15 for all the pellets. Hence 
we can use Eq.(5) from Mennell a et alJ (1199 8) to calculate g(e, 
f). We have adopted the same value as lMennella et al.l (119981) for 
the dielectric function of the host mediuni (PE) and for the grain 
material. However, adopting lBoudet et al.l ([2005) values does not 
changes the value of g(e, f) by more than 2-3%. We find that g(e, 
f) is usually greater than 0.9 except for pellets containing 50 and 
100 mg of samples for which it is equal to 0.78. The transmit- 
tance corrected for reflexion loss, Tcorr, is calculated using the 
formula: T^orr = Tnieasured/(1 " R)^ with R, the reflectance of the 
pellets, expressed as R = (1 - neff)^/(l + ngfj)^. The effective re- 
fractive index of the pellet, n^ff, is related to the average dielec- 
tric function of the pellet, feff, by the formula nej - ^Re{eeff). 
The correction factor for reflexion loss, (1-R)^, is larger than 0.9 
for all our pellets. 

Once these corrections have been applied to the different 
spectra corresponding to different amount of samples and dif- 
ferent experimental configurations, they are merged into a sin- 
gle spectrum by taking advantage of the overlapping regions of 
the different spectra. The spectra were merged with respect to 
short wavelengths measurements. Typically the absolute opac- 
ity of the different spectra to be merged agree within less than 
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Fig. 2. Mass absorption coefficient of tfie syntfiesised samples in tfie mid infrared range. Panel (a): amorphous samples with olivine 
stoichiometry, Fl, F2, F3, panel (b) crystallised Fl, F2 and F3 samples, panel (c) amorphous samples with enstatite and diopside 
stoichiometry, E and D, respectively, panel (d) crystallised E and D samples. In panel(d) the absorption-like feature marked by an 
asterisk is an experimental artefact. 



10% at short wavelengths and within ~ 10-30% in the 250-300 
fim range, i.e. between experimental configurations using the 23 
and 50 jum Mylar beamsplitters. Thus we estimate that the un- 
certainty on the absolute value of the MAC, k, is of the order 
of 10 to 30%. However, it is important to note that, the rela- 
tive uncertainty on the MAC of a sample measured at different 
temperatures within an experimental run is much lower, of the 
order of 5% at most. This uncertainty is mainly due to the over- 
all stability of the spectrometers which is critical at very long 
wavelength. The spectrometer stability was evaluated by repro- 
ducibility measurements and by comparing overlapping regions 
of spectra taken with different configurations and also of spectra 
taken with different spectrometers (the one from the ESPOIRS 
setup and the one from the AILES beam line at the SOLEIL 
synchrotron). The measurements at very long wavelengths were 
performed on the ESPOIRS experimental set-up up to 800 yum 
and were merged with measurements performed on the AILES 
beam line at the SOLEIL synchrotron in the range 500-1500 
fim. Above 800 fim the measurements were performed solely 
with the AILES beam line setup. However, thanks to the high 
brightness and stability of the FIR synchrotron radiation and to 
the good performance of the spectrometer, the relative uncer- 
tainty on the transmission does not exceed 5 %. This was tested 



from measurements performed during runs separated by several 
months. It is co mpar able to the trans mission accur a cy sta ted by 
Me nnella et al.i (Il998h . 5-10 %, and lBoudetet"an (l2005h . 2%. 
The absolute uncertainty on k is expressed as follow: Sk - S/M 
xST/T, where dT/T is the relative uncertainty on the measured 
transmittance, which is estimated to 2.5 % and 5 % for spectra 
taken with the 23 jum and with the 50 or 125 fim beamsplitters, 
respectively. 

4. Results 

4.1. Mid Infrared spectra 

We report in Fig. |2] panel a and c, the MAC spectra of all the 
amorphous silicate samples in the 7-27 yum range. We observe 
the two characteristic bands around 10 fj.m and 20 jum which are 
attributed respectively to the Si-O stretching and O-Si-O bending 
vibrational modes. For samples Fl, F2, E and D, the two bands 
are broad and featureless confirming the amorphous nature of 
these samples. However, unlike the MIR spectra of the Fl and 
F2 samples, the MIR spectrum of F3 clearly shows vibrational 
bands from crystalline material (at ~ 11, 16, 19 jjm) superim- 
posed onto the broad and featurelesslO and 20 yum bands. Thus, 
the F3 sample is a mixture of amorphous silicate with forsterite 
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Fig. 3. Mass absorption coefficient of amorphous olivine samples at different temperatures in the 100-1000(1200) fim range. Panel 
(a): sample Fl: Mg2.3Si04, panel (b): sample F2: Mg2.8Si04, panel (c) sample F3: Mg2.05SiO4. 



whereas all the other silicates are completely amorphous. From 
the ^^Si NMR spectrum of this sample, we estimate the amor- 
phous crystalUne fraction to be ~ 14% of the sample mass (see 
SectO- 

The bands from forsterite dominate the spectra of the Fl, 
F2 and F3 samples annealed at 1100 °C (Fig. |2l panel b). The 
spectra of the three samples Flc, F2c and F3c are very similar 
and representative of the forsterite spectrum. The position and 
strength of the main spectral features of sample Flc, F2c and 
F3c are i n agreement with t he spectra of natural or synthesised 
forsterite (iFabian et al.ll2001 ). This is true even for the Flc sam- 
ple which contains a significant amount of enstatite (see Table|2l). 
The spectra of the E and D sample annealed at 1 100 °C, Ec and 
Dc (Fig. 121 panel d) are dominated by the characteristi c vibra- 
tional bands of respectively enstatite and diopside (Koike et alJ 
l2000h . 

4.2. Far Infrared and temperature dependent spectra 

This section describes the main results obtained from the 
spectroscopic measurements performed in the FIR/submm 
rang43- As exposed in Sect. [T] and discussed in more detail in 



' The mass absorption coefficients of the studied samples are avail- 
able upon request to Karine Demyk (demyk@cesr.fr). 



Sect. |6] simple semi-classical models allow us to express the 
grain MAC in the FIR/submm spectral range with the following 
formula (same as Eq.|2]): 

K^K(Ao)d)-f' (5) 

AO 

The parameter /3, which is the slope of the MAC, k, in 
logarithmic scale, in the FIR/submm range, is similar to the 
emissivity spectral index used in astronomical studies. These 
models predict a single value for yS in the FIR/submm domain 
that is independent of the temperature. Although /? is not a usual 
parameter to describe experimental data, we use it in the present 
study for comparison with the astronomical emissivity spectral 
index. 

Samples Fl, F2 and F3 (Mgz.Qs isSiOi): 
Fig. [3] shows the evolution of the MAC, k, of the amorphous 
samples Fl, F2 and of the partially crystalline sample F3 as a 
function of A in logarithmic scale and for different grain tem- 
peratures. Commonly for these three samples having an olivine 
stoichiometry, we observe that ( i) the MAC diminishes when the 
grain temperature decreases, (ii) this diminution is more pro- 
nounced at long wavelengths and (Hi) the shape of the MAC 
varies with the wavelength and its slope gets steeper at long 
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Fig. 4. Mass absorption coefficient of amorphous samples with pyroxene stoichiometry at different temperatures in the 100- 
1000(1500) fim range. Panel (a): E sample, Mgo.gsSiOs, panel (b): D sample, Cao.ggMgo.gSiaOg. 



wavelengths. For instance, the MAC at 1 mm, decreases by a 
factor of 2.1, 3.5 and 3.8 from 300 to 10 K, for the F2, F3 and 
Fl samples, respectively, whereas the decrease is negligible at 
110 fj.m. We note that, for the 3 samples, the MAC does not vary 
for temperatures smaller than 30 K, the MAC measured at 10 
and 30 K being almost indistinguishable. For all temperatures, 
the spectra show an inflexion in the absorption curves that is 
more or less pronounced, and that extend on a variable spectral 
domain depending on the considered sample. This break occurs 
around 550-700 yum depending on the sample but also depend- 
ing on the grain temperature (it occurs at shorter wavelength at 
high temperature). It is therefore impossible to use Eq. |5]to fit 
the MAC of these samples since a single spectral index cannot 
be defined over the entire FIR/submm range. However, on both 
sides of this inflexion point a linear behaviour is observed and we 
can thus define two spectral regions in which we can calculate 
j6, the slope of the curve (Table[3]l. For the Fl sample, at 300 K, 
the spectral index calculated in the 130-670 fim range is = 1.6 
whereas in the 710-1200 /urn range, /3 - 2.2. At 10 K, /3 is equal 
to 2.1 in the 130-690 yum range whereas it is 3.6 in the 690-1200 
fj.m range. For the F2 sample at 300 K, in the range 170-770 fim, 
P = 1.8 whereas jS = 2.5 in the range 770-1000 yum. At 10 K, in 
the same wavelength ranges, yS is respectively equal to 2.1 and 
3.2. For the F3 sample at 300 K, the spectral index calculated in 
the range 140-550 yum is equal to 1.9 whereas in the 550-1200 
yum range it is equal to 2.5. At 10 K, /? is respectively equal to 1 .9 
and 4.5 in the wavelength ranges 150-650 //m and 650-1200yL(m. 

Even though all the three samples present a similar fre- 
quency and temperature dependence, some differences exist 
between their MAC. In the case of the Fl and F2 samples, we 
observe that the MAC decreases with decreasing temperature 
on the whole wavelength range whereas for the F3 sample, the 
variation of /<• as a function of temperature is present only for 
A > 550-650 //m. This may be due to presence of crystalline 
material, forsterite, in the F3 sample. Interestingly, despite the 
crystallised fraction present in the F3 sample (14% in mass, see 
Sects. 12.51 and 12.41 ) and that is clearly observable in the MIR 
spectrum of this sample, its FIR optical behaviour is dominated 
by the 86% amorphous fraction, as shown by the variation of 
its MAC with the temperature. The minor differences found for 
the values of fi, for a given wavelength range and at a given 
temperature, between the remaining 2 samples Fl and F2 may 



be due to differences in composition (Mg excess), homogeneity 
and microstructure induced by the sol-gel methods used for the 
synthesis (see Sect.|2|. 

Sample E (Mgo ^sSiOi): 

For the E sample, we observe that (/) the mass absorption 
coefficient decreases when the grain temperature decreases, 
the diminution of the MAC with the temperature is more 
pronounced at A > 400 fj.m, and the shape of the mass 
absorption coefficient varies with the wavelength and its slope 
gets flatter at long wavelengths (Fig. H) panel a). Indeed, we 
note that at 200 fim, the mass absorption coefficient measured at 
10 K is lowered by a factor of 2.2 with respect to the value at 
300 K. This factor is equal to 3.9 at 400 yum, to 6.8 at 800 fim 
and remains constant up to 1500 fim. In contrast to the F 
samples, the MAC of the E sample varies significantly between 
10 and 30 K. As for the F samples, the shape of the MAC 
clearly does not follow the asymptotic behaviour predicted by 
the semi-classical models: a single /3 value cannot reproduce 
the experimental spectra. Unlike the F samples for which the 
MAC gets steeper, the MAC of the E sample is flattened at long 
wavelength. It is not possible for this sample, to define only two 
regions characterised by a linear behaviour with the wavelength. 
Even though it is more difficult to define separated spectral 
regions than for the F samples, we identify three regions in the 
spectra for which we can define and calculate a spectral index. 
Each region is delineated by an inflection in the absorption 
curves at 420-590 jum and at 790-800 jum, depending on the 
grain temperature. The spectral index, /3, calculated in these 
3 wavelength ranges is indicated in Table [3] We note that at 
a given grain temperature, /3 decreases when A increases. The 
spectral index derived from the MAC measured at 300 K is y6 = 
1.7 in the 150-420yum range, it decreases to 1.3 in the 420-800 
yum range and to 0.9 in the 800-1500 fj.m range. We observe 
the same flattening for the spectral index derived from the 
MAC measured at 10 K: yS = 2.5, 1.7 and 0.9 in the wavelength 
ranges: 150-590, 590-790 and 790-1500 //m, respectively. At a 
given wavelength, /3 increases when the temperature decreases 
only for the wavelength range < 800 yum. Above 800 /vm, /3 is 
constant with the temperature. 

The MAC of the E sample i s in qualitativ e agre ement with the 
MgSi03 sample studied by iBoudet et alJ (l2005i) even though 



8 



A. Coupeaud et al.: Low temperature FIR and submm mass absorption coefficient of interstellar silicate dust analogues 

Table 3. Value of the spectral index p derived on different spectral ranges from the experimental mass absorption coefficient spectra 
for the studied amorphous samples. 



Sample 


spectral domain (//m)^^* 


P (10 K) 




p (lUU K) 


Q /inn 
p (2UU Kj 


p \iW JS-) 


Fl (Mg2.3Si04) 


130 - 690/710 


2.1 


2.1 


2.0 


1.8 


1.6 




690/710- 1200 


3.6 


3.8 


3.4 


3.2 


2.2 


bZ (Mg2.8alU4) 


1 /U - / /U/oUU 


O 1 

Z.l 


L.L 


z.l 


Z.U 


1 Q 
i.O 




770/800 - 1000 


3.2 


3.2 


2.9 


3.0 


2.5 


F3 (Mg2.05SiO4) 


150 - 550/650 


1.9 


1.9 


1.9 


1.9 


1.9 




550/650 - 1200 


4.5 


4.5 


3.4 


3.0 


2.5 


E(Mgo.95Si03) 


150 - 420/590 


2.5 


2.2 


2.1 


2.0 


1.7 




420/550 - 800 


1.7 


1.5 


1.5 


1.4 


1.3 




800 - 1500 


0.9 


0.9 


0.9 


0.9 


0.9 


D (Cao.98Mgo.9Si206) 


150 - 450 


2.0 


2.1 


2.0 


1.9 


1.8 




450 - 650 


1.4 


1.4 


1.3 


1.3 


1.3 




650 - 1000 


0.8 


0.7 


0.6 


0.6 


0.9 



Notes. 

The number of spectral domains on which a spectral index, p, may be defined depends on the samples. In addition, for a given sample, the 
spectral domains may vary with the grain temperature. When two values of the wavelength separating the spectral domains are given, the first one 
corresponds to the highest temperature (300 K) and the second one to the lowest temperature (10 K), see text. 



some differences exist that may be attributed to differences 
in the samples arising from different synthesis methods. In 
both experiments, the MAC is characterised by a flattening 
at long wavelength. However, the comparison between the 
two measurements is not straightforward. Indeed, the values 
iBoudet et al.l (l2005h derived for p are not comparable with ours 
because they are calculated in two different spectral ranges 
(100-200 ^xm and 500-lOOOyum) that are not adapted to the 
sample E presented in this study. This illustrates the problem of 
the relevance of the use of the spectral index j3 (see Sect l6.2l i. 

Sample D (Caft tf^MgQ ^SiiOf,): 

The variations of the MAC of the D sample with the wavelength 
and with the grain temperature (Fig. [H panel b) is qualitatively 
similar to what is observed for the E sample: (i) decrease of 
the mass absorption coefficient with grain temperature that ( ii) 
is more pronounced at long wavelength, and (Hi) flattening of 
the MAC at long wavelength. However, quantitatively, these 
variations are much less pronounced in the case of the D 
sample. For the amorphous D sample, we observe the decrease 
of K with decreasing temperature but this diminution is almost 
constant above ~ 250 fim. As for the other F and E samples, the 
MAC cannot be described by a single spectral index p. For all 
temperatures, the spectra are characterised by an inflection point 
in the absorption curves at ~ 450 fj.m and another at ~ 650 yum. 
Like for the E sample, these inflection points correspond to 
successive flattening of the MAC. Indeed, between 150 and 
450 fim, we calculate at 300 K a spectral index with a value 
of 1.8 that decreases to 1.3 between 400 and 600 fim and then 
to 0.9 in the 650-1000 fim. Finally, we observe at a given 
wavelength, that p is constant with the temperature (Table |3]l. 

Crystallised Flc, F2c, F3c, Ec and Dc samples: 

To investigate the hypothesis that the amorphous state could be 
responsible for the variation of the grain MAC as a function of 
temperature and wavelength, part of each sample were heated 
for 24 hours at 1100°C to crystallise them. At this temperature 
we know from the MIR (Fig |2]i and the NMR and XRD analysis 
that our samples are crystalline. We report the evolution of k as 



a function of A in logarithmic scale at different temperatures for 
the crystallised samples Flc and F3c (olivine stoichiometry) and 
for the crystallised samples Ec and Dc (pyroxene stoichiometry) 
(Fig. |5]l. The shape of the mass absorption coefficient deviates 
from the asymptotic form characterised by a single value for p. 
This may be due to the mixture of the different crystalline phases 
that are present in the annealed samples. During the process of 
crystallisation, separate phases have been formed. The F sam- 
ples have crystallised in 2 or 3 phases depending on the sam- 
ple (see Sect 12.41 and Tabled: forsterite (Mg2Si04), periclase 
(MgO) and ortho-enstatite (MgSiOs), the E sample in 3 phases: 
ortho-enstatite, forsterite, and clino-enstatite and the D sample 
in 3 phases: diopside, akermatite and merwinite (see Table |2]l. 
However it is important to note that, in contrast to the amor- 
phous samples, the MAC of the crystallised samples does not 
vary with the temperature. The MAC value stays constant for 
all grain temperatures. In addition its spectral shape remains un- 
changed for all grain temperatures, i.e. the spectral index does 
not vary with the grain temperature. The optical properties of 
amorphous and crystallised samples thus behave very differently 
in the FIR/submm range as a function of the temperature. 



5. Discussion 

The study of the dust analogue MAC in the FIR/submm and 
of its temperature dependence reveals differences between the 
various materials (samples having an olivine stoichiometry: 
Mg2.()5-2.8Si04 versus samples with a pyroxene-stoichiometry: 
Mgo.gsSiOs and Ca().98Mgo.9Si206), and at the same time a sim- 
ilar behaviour. The common behaviour can be characterised as 
follow: 

1 . All amorphous samples present a MAC that increases when 
the temperature rises, and this variation vanishes when the 
samples are crystallised by annealing. This implies that 
thermally-activated processes causing an additional absorp- 
tion exist only in the amorphous state. When the tempera- 
ture increases, the thermally activated processes increase the 
MAC, first in the long wavelength range (A > 500-600 yum) 
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Fig. 5. Mass absorption coefficient of crystallised olivine and pyroxene samples at different temperatures in the 100-1000(1200) fim 
range. Panels (a-b): Mg2Si04 samples, panel (c) MgSiOs and panel (d) CaMgSiaOe. 



at low temperatures and then progressively over the whole 
FlR/mm range at higher temperatures. 

2. At the lowest experimental temperatures, when these 
thermally-activated processes do not add a significant 
absorption, the remaining temperature-independent MAC 
presents a common behaviour in all our amorphous samples: 
in the short wavelength range (100 < A < 500-700 jum), the 
local spectral index, has a value of the order of 2, whereas at 
longer wavelength, the value of /3 varies after a more or less 
abrupt wavelength transition (Table O. 

The differences in behaviour of the MAC of the amorphous sam- 
ples of different chemical compositions (olivine and pyroxene 
composition) with the wavelength and temperature are the fol- 
lowing: 

3. The Fl, F2, and F3 samples (olivine stoichiometry) present 
a MAC characterised by a spectral index around 2, 
temperature-independent in the 120-650 yum range, followed 
by a steep transition towards spectral index values greater 
than 3 at /I > 700 fim. Despite the differences in composition 
due to different synthesis recipes and despite the presence of 
a small amount of crystalline material in one of the sample, 
the behaviour of the 3 samples with an olivine composition 
samples are similar. 



4. The behaviour of the 2 samples with a pyroxene composi- 
tion, E and D, is different from the one of the samples with 
an olivine composition. The spectral index value around 2 in 
the short wavelength range evolves with a smooth transition 
towards lower values, down to 1 for A > 800 fj.m. However, 
despite this similar behaviour, the variation of the MAC 
with the temperature is much higher in the E sample than in 
the D sample. 

From point (1) it is clear that the variation of the MAC with 
the grain temperature is related to the disorder of the mate- 
rial, i.e. to its amorphous nature. This is clearly shown by the 
disappearance of these variations when the materials are crys- 
tallised, i.e. well-ordered. This is true for all samples: the Fl, 
F2, and F3 samples with an olivine stoichiometry, the E and 
D samples with a pyroxene stoichiometry as well as for siUca 
samples (Coupeaud et al. in preparation). A similar behaviour 
was also observed by Mennella et al. ( 1998) for iron-rich amor- 
phous and crystalline samples with an olivine stoichiometry 
(Mgi i2Fei 88Si04). Hence, the MAC variation with the grain 
temperature is a general characteristic of amorphous materi- 
als, regardless of their composition or detailed micro-structure. 
However, the amplitude and wavelength dependence of this 
variation is different for the various samples (olivine- versus 
pyroxene- stoichiometric samples but also within the pyroxene- 
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or olivine stoichiometric samples). The origins of these differ- 
ences are likely to be linked with the structural state of the amor- 
phous samples which differs by their compositions and prepara- 
tion methods. The detailed description of an amorphous phase 
is difficult as well as to identify measurable physical parame- 
ters that may be responsible for the observed MAC variations. 
Crystals are solely described by their composition and lattice 
periodicity but this is not the case for amorphous solids that are 
characterised by disorder of chemical nature (at a given position 
into the solid the composition of the atom changes) and of geo- 
metrical nature (the bond lengths and angles vary). Even though 
they may have broad structural order at the short (nearest neigh- 
bours) and intermediate (second or higher coordination sphere) 
length scales, amorphous solids are characterised by the lack of 
long-range periodicity. Consequently, in amorphous solids there 
is flexibility in the way that atoms can be organised in a variety 
of local arrangements, thus resulting in a diversity of structures 
at short and medium length scales. Composition is certainly a 
parameter that may be used to describe amorphous solids since 
it has a strong influence on the atomic arrangement at short- and 
intermediate-range order but it is not the only one since amor- 
phous materials having the same composition may have differ- 
ent microscopic structures as this is the case for glass and solgel 
materials. 

If one considers only the spectral shape of the MAC of the stud- 
ied samples at a given grain temperature, the experiments pre- 
sented in Sect. H show that the FIR/submm MAC of the Fl, F2 
and F3 samples and of the E and D samples are different (points 3 
and 4). Interestingly, experiments performed on silica samples, 
Si02, have revealed a spectral shape of the MAC of silica dif- 
ferent from the olivine and pyroxene samples (Coupeaud et al, 
in preparation). Hence, this seems to indicate that the spectral 
shape of the MAC of amorphous silicate samples depends, at 
first order, on the cation/Si ratio of the material, which is equal 
to for the siHca samples, (Mg2++Ca2+)/Si = 0.95 and 0.94 
for the samples with a pyroxene stoichiometry, E and D sam- 
ples, respectively, and Mg^^/Si = 2.05, 2.3 and 2,8 for the F3, 
Fl and F2, samples with an olivine stoichiometry, respectively. 
The cation/Si ratio seems to be a better parameter than the com- 
position, since the samples E and D, have different composi- 
tion but their MAC has qualitatively the same spectral shape. 
This is also true for the 3 Mg-rich samples which MAC has 
also qualitatively the same spectral shape as the one of Fe-rich 
oli vine-stoichiometry FA YA sample (Mgo.igFei g2Si04) studied 
bv lMennellaetalJ ( ll998l) . However, because iron is an important 
component of interstellar dust it would be interesting to study 
more Fe-rich samples (with iron within the silicate network for 
different stoichiometrics from enstatite to olivine, with oxide in- 
clusions) to investigate whether its presence within the grains 
may induce or not a different spectral behaviour of the MAC. 
However, for a given value of the cation/Si ratio, the ampli- 
tude of the MAC variations changes from one sample to an- 
other as this is the case for the E and D samples. These 2 sam- 
ples have different compositions, with the D sample contain- 
ing Ca^"^ cations in addition to Mg^^ cations. The presence of 
Ca^^ cations, which are linked to 6 oxygen atoms whereas Mg^^ 
cations are linked to 4-5 oxygen atoms, modifies the chemical 
and geometrical microstructure of the D sample compared with 
the E sample and it certainly participates to explain the MAC 
variations. As far as the samples with an olivine stoichiometry 
are concerned, the differences between the 3 samples may also 
be explained by differences in their stoechiometry (the Mg^^/Si 
ratio varies from 2.0 to 2.8) and thus in their chemical micro- 
scopic structures resulting from the use of different recipes for 



the solgel synthesis. Indeed the XRD analysis show that the crys- 
tallised samples Flc, F2c and F3c are a mixture of different crys- 
talline phases and it is likely that the amorphous samples contain 
precursors of these phases such as for example chains of Si04 
tetrahedra for the Fl sample (precursor of enstatite), more or 
less MgO inclusions in samples Fl and F2 and partly crystallised 
material in the F3 sample. 

The observed variations of the MAC may also be understood in 
term of degree of connectivity of the Si04 tetrahedra in the ma- 
terial (obtained from the NMR analysis). The samples having 
an olivine stoichiometry contain a majority of isolated tetrahe- 
dra (Q'^ species) and dimers of Si04 (Q' species) whereas the 
samples having a pyroxene stoichiometry contain longer chains 
(Q^ species) and 3-dimension structures (Q"* and Q"* species). 
The sample having a stoichiometry of diopside, for which the 
MAC variation is smaller than for the sample having an enstatite 
stoichiometry, contains much more and Q"* species than the 
sample having an enstatite stoichiometry. The link between the 
connectivity of the material and the optical behaviour at low tem- 
perature should be investigated further via more sensitive and 
detailed NMR studies. 

Our samples, because of the synthesis method, contain OH 
groups and water that may influence the MAC variation of 
the silicates samples in a similar way a s for silica samples 
(Agladze et al. 1996: Boud et et al.ll2005h . iBoudet et alJ (l2005h 
proposed that OH groups could be involved into tunnelUng pro- 
cesses in two-level systems (TLS) explaining partly the exper- 
imental MAC variations in Si02 samples. From the OH vibra- 
tional bands at ~ 3 fim, we estimate that all samples have roughly 
the same amount of water and OH groups. Hence, even though 
water and OH groups may participate to the temperature depen- 
dance of the MAC, they are certainly not responsible for the ob- 
served spectral shape of the MAC since the wavelength depen- 
dence of the MAC of the F samples, with an olivine stoichiome- 
try, is different from that of the E and D samples, with a pyroxene 
stoichiometry. As said above it is probably related to the cations 
(Mg^^ and Ca^^) and to their abundance with respect to silicon. 
It is interesting to note that the different spectral behaviours in 
the FIR range of the F samples compare to the E and D sam- 
ples are accompanied by different absorption properties in the 
8-12 fim range, i.e. in the Si-O stretching mode of the samples. 
The stretching modes of the E and D samples are very simi- 
lar and peaks toward shorter wavelengths compared with the 
F sample stretching mode (Fig. |2]l. Hence, the MIR spectrum 
could be used to predict the FIR behaviour of the dust analogues. 
Unfortunately this may not hold for the astronomical silicates 
for which the MIR spectrum most certainly reflect a mixture of 
amorphous silicates with various compositions. 
The common behaviour of all samples (points 1 and 2) comes 
from the general properties characteristic of the amorphous 
state. It may be understood in the framework of the TLS model 
developed within our group (Menv et al. 2007). The TLS model 
decompose the MAC in two components, one that depends 
on the grain temperature and the other that is temperature- 
independent. It predicts that the temperature-independent MAC 
in the FIR/mm spectral range is characterised by a local spectral 
index varying from y6 = 2 at short wavelengths to = 4 at 
long wavelengths. The transition between these two values 
usually occurs in the 500-1000 jjm range, depending on a 
characteristic length that describes the charge disorder in the 
amorphous material. This temperature-independent behaviour 
is expected to be observed at low temperature, typically, T ~ 
10 K. When the temperature increases, a thermally activated 
additional opacity appears, characterised by a lower value of 
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the local spectral index that may be as low as 1 . The amplitude 
of the variation of this additional opacity with the temperature 
depends on the amount of disorder in the amorphous structure 
and this disorder-induced opacity disappears when the sample 
is crystallised. When the temperature is high enough, this 
thermally activated opacity should dominate and determine the 
spectral shape of the opacity of a given amorphous material. In 
the TLS model, the thermally activated opacity is due to pres- 
ence within the amorphous structure of closed rot ational-type 
config urations of atoms or groups of atoms (see iMenv et aU 
(l2007h for a more detailed description). These may be local 
defects or impurities. It has been observed for example in 
various silica samples in which the presence of OH-groups is 
respo n sible for the therm ally activated opacity (lAgladze et alj 
[19961 lBoudetetalJl200l . It may also be groups of atoms 
participating to the material structure itself. The fact that in 
our experiments amorphous samples with olivine and pyroxene 
stoichiometry have different opacity variations suggests that 
network modifiers such as Mg^^ and Ca^^ ions do not have well 
defined positions in the amorphous structure, allowing such 
closed configurations to exist. 



6. Astrophysical Implications 

6.1. Relevance of the dust analogues studies 

Experimental data, even on a limited number of samples, are 
always relevant because they represent a real behaviour which 
cannot be ignored. It consists in the present case in the fact 
that the opacity of amorphous silicate grains varies with tem- 
perature, and cannot be described at any temperature with a 
single spectral index over the whole FIR/mm range. These 
experimental results are compatible with the anti-correlation 
between the emission spectral index and the temperature re- 
vealed by astrophysical observations from experiments such as 
PRONAOS, Arche o ps, Herschel or Pla nck (Dupac et al. 20031 
Desert et al.' 12008'; R odon et al.1 120101: [Anderson et al.i i20ia 
Veneziani et al. 2010) and with the submilli metre excess ob- 
served by FIRAS/WMAP ( Paradis et all2009h . However, at the 
same time experimental data on interstellar dust analogues must 
be carefully considered. The comparison of the experimental 
data with the observations is not straightforward since the astro- 
nomical FIR/submm data average the emission from dust grains 
that may be coming from diff'erent regions of the ISM and/or 
that may have different sizes, shapes, chemical compositions and 
temperatures. 

The studied silicates grains have an uncharacterised size distri- 
bution with typical grain sizes of a few tens to a fews hundreds 
of nm. Interstellar silicate grains, that may have a grain size dis- 
tribution different from the one of dust analogues, have a typi- 
cal maximum size of a few hundreds of nm, much smaller than 
the wavelength in the far infrared. Experiments on silica spher- 
ical grains with diameters ranging from 7 to 500 nm show that 
the measured grain opacity does not depend on the grain size 
(Coupeaud et al. in preparation). In addition, the absorption part 
of the refractive index being low, the Rayleigh criteria is ful- 
filled, and thus the grain size distribution should not influence the 
resulting dust opacity spectrum. Even though the shape distribu- 
tion of interstellar grains is unknown, there is no reason for them 
to be spherical. The polarisation observed in the 9.7 fim band 
implies that pa rt of the grains are aligned and thus not spherical 
(IWhittetll2003K Interstellar grains could have all possible forms, 
although they may be rounded by erosion induced by sputtering. 



The opacities presented in this study are those of a population of 
submicronic grains having a natural shape distribution resulting 
from the synthesis method used. They should thus be relatively 
close to what can be found in the ISM. Interstellar grains may 
also coagulate to form dust aggregates, in which case the shape, 
the fluffyness and the properties of the aggregates also need to 
be taken into account. Experimental studies of dust aggregates 
are difficult to perform but modelling show that dust aggregation 
increases the dust emissivity without changing the shape of the 
emissivity as a function of wavelength (see Sec. 16.31 ). The stud- 
ied samples thus represent relevant analogues of interstellar dust 
in terms of composition but also in terms of size and shape. 

6.2. Classical dust emission models and relevance of the 
spectral index 

Because no specific mechanism for absorption and emission was 
expected in the FIR/mm range, the opacity in that domain was 
considered to be dominated by the long wavelength wings of 
the observed mid-infrared vibrational bands. The temperature- 
independent Lorentz model for bound oscillators is then well 
appropriate. In the FIR/mm range, far from the MIR resonances, 
it leads to the simple asymptotic behaviour for the dust opac- 
ity indicated Eq. |5]and for the dust emission for optically thin 
medium (Eq.[T]) to the following expression: 

I(^) = B(A, T) MM)W^or'' d-2 (6) 

where p - 2. This model fully describes the emission with 
only two dust parameters: the opacity k(Aq) at a wavelength, 
Aq, where the asymptotic behaviour of the Lorentz model 
is assumed to be fulfilled, and the dust temperature, T, that 
depends on its radiative environment. The Lorentz model 
applies for crystals, but it is also used to describe the FIR/mm 
properties of amorphous solids. The mid-infrared bands, broad- 
ened by the dispersion of angles in the amorphous structure, 
are usually modelled by the sum of closed Lorentz oscillators 
(up to 33 osc il lators are considered in the astrosilicates from 
iDraine & Led (1984)). This modelling does not change the 
asymptotic behaviour in the long wavelength range which is the 
same as for crystals. The dust emission may also be explained 
in terms of phonons (iHennin^ & Mutschkelll99'7l) . Because of 
the breakdown of the momentum selection rule in amorphous 
materials, each phonon can be involved in the temperature- 
independent (or weakly dependent) one-phonon absorption 
process. This disorder-induced one phonon process gives rise 
in the FIR/submm to a continuous absorption that completely 
dominates over multi-phonon absorption processes. As a con- 
sequence the dependence of the opacity with the wavelength is 
expected to follow the one of the phonon density of state g{A), 
which follows a simple power law dependence in the Debye 
approximation: g(A) oc A'^ in 3-dimensional structures, and 
g(A) oc in 2-dimensional structures (Mennella et al. 199^. 
These two models thus predict that the FIR/mm dust-emission 
follows at first order the law given by Eq.[6|in which the opacity 
is temperature-independent and where yS equals 2, or may 
vary between 1 and 2. This variation of the spectral index is 
usually understood in terms of different grain properties: some 
atomic structures or surface-to-volume ratio which can induce a 
2-dimensional behaviour, and/or a distribution of temperatures 
along the line of sight to justify a decrease of the spectral index 
to values smaller than 2 (ISeki & Yamamoto[ll980l) . 

Within the framework of these models, astronomical ob- 
servations are thus interpreted with a single, temperature- 
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Fig. 6. Dust emission simulated using the modified blackbody 
law (see Eq|6]l for lOK interstellar grains. The solid lines are 
calculated using the experimental lOK opacity of the amor- 
phous E (MgSiOs) sample (black line) and the amorphous Fl 
(Mg2.3Si04) sample (grey line). The dust emission assuming 
a single spectral index of /3-2.0 (dotted line), ^6=3.4 (dotted- 
dashed line) and j6= 1 . 1 (dashed) are also represented. The curves 
have been normaUzed at 600 yum. 



Fig. 7. Dust emission simulated using the modified blackbody 
law (see Eq|6]l for a mixture of amorphous grains with enstatite 
and olivine stoichiometry, for two grains temperatures: 10 and 
30 K. The solid grey lines are calculated using a single spectral 
index, j6=2. The amorphous E sample with an pyroxene stoi- 
chiometry and the Fl sample with an olivine stoichiometry 
have been used for the modelling with the following respective abun- 
dances: 0.2:0.8 (dashed line), 0.5:0.5 (dotted line) and 0.8:0.2 (dotted- 
dashed line). The curves have been normalized to unity. 



independent, spectral index value over the whole FIR/mm range. 
As explained above, such an interpretation relies on the assump- 
tion that the knowledge of the dust opacity at a given FIR- 
wavelength Aq (deduced from laboratory measurements at any 
temperature) and the dust spectral index value allows us to de- 
duce the emission of the dust at any FlR/mm wavelength and 
at any temperature. However it is clear from the experimental 
results above that in the amorphous state, the grain opacity, a:, 
varies with the temperature and the wavelength. Therefore it is 
not possible to describe the dust opacity with a single spectral 
index. It may then be more appropriate to use the local spectral 
index, /5iA), define as the slope of the opacity at any wavelength, 
to take into account the opacity variations with the wavelength. 
This has strong implications that question the relevance and the 
physical meaning of the emissivity spectral index derived from 
the observations. First of all it means that the use of a single 
spectral index should not be the rule when modelling the dust 
emission. Indeed, since the experimental spectral index value of 
relevant interstellar dust analogues varies with the wavelength, 
the emissivity spectral index deduced from astronomical data 
may also vary with the wavelength. As a consequence, obser- 
vations of an astronomical object at different wavelengths in the 
FIR/mm range can lead to different values of the emissivity spec- 
tral index (and slight variations in the dust temperature). It is 
inter esting to note h ere that the FIR opacity of the astrosilicate 
from lDraine & Led (11984) was modified to account for the high 
galactic latitude dust emission measured by FIRAS and that, t he 
FIR opacity of the new astrosilicates fromlLi (S:Drainel(l2001l) IS 
wavelength-dependent, the spectral index being equal to 2 below 
800 jjm and to 1 .6 above. Secondly, it means that the knowledge 
of the opacity at a single FIR wavelength and a single spectral 
index, either expected from simple theory (J3=2) or measured in 
a given limited wavelength range, is not sufficient to derive ac- 
curately the emission of the dust on the whole FIR/mm domain. 
These two points are well-illustrated by Table[3]and Fig.|6]that 
shows the simulated emission of two amorphous grain popula- 



tions of mean composition close to MgSi03 and Mg2Si04. For 
A < 600 fim, the derived spectral index is/3 ~2 whereas it is ~ 4 
at A > 600 yum for the grains with an olivine stoichiometry and, 
respectively,/? ~ 2.5 and ~ 1.1 for the grains with an enstatite 
stoichoimetry. 

However, the fact that the dust emission in certain astronomical 
environments are well-modelled with a single t emperat ure and 
emissi vity spectral index value ( Planck Collaboration et al.l 
1201 lal) is not in contradiction with the experimental results 
presented in this study. In this particular case, a population of 
amorphous grains with olivine and pyroxene stoichiometry may 
produce a resulting emissivity spectral index constant over the 
whole spectral range (Fig. |7]i. Actually, such a case illustrates 
well the difference between the experimental spectral index and 
the observational emissivity spectral index that must be care- 
fully used to describe the dust as it averages contributions from 
various dust composition. It also illustrates that the variation of 
the dust spectral index towards various lines of sight may trace 
changes in the dust composition. 



6.3. Dust mass and opacity 

These experimental results have also consequences on the use 
of methods based on flux ratios to determine the dust temper- 
ature and spectral index that ar e comm only used in absence of 
spectroscopic data. lShettv et al.l (l2009h investigated the effect of 
line of sight temperature variation and noise on the accuracy of 
the derived parameters such as the dust temperature and spectral 
index or the dust mass. However, another important source of 
error comes from the fact that the flux ratios method is based on 
a single value for the spectral index and from uncertainties on 
the MAC itself. From the experimental results presented here it 
is also possible to investigate this effect. Flux ratio are propor- 
tional to the ratio of the blackbody emissions at the dust temper- 
ature and to the ratio of the opacities, at the two wavelengths. 
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The error may be important if the fluxes are observed at wave- 
lengths that are very different and situated on both sides of ~ 
600-800 i-im, where the local spectral index of dust analogues 
changes. As an example the ratio of the experimental MAC at 
450 and 1200 /im for the E sample (pyroxene stoichiometry) is 
'<'450 / *'i200 ~ 4.8 and 3.7 at lOK and 30 K, respectively, whereas, 
using the asymptotic expression for the dust MAC (Eq|5]l, /C45o / 
Kiioo ~ 7. 1 for /3-2 and K450 / /<-i2()o ~ 11 .6 for j6 = 2.5 (value that 
reproduces the MAC of the E sample at /I < 600 fim). Hence, the 
flux ratio is wrong by a factor of 1 .5-2 if one adopt /3=2 and by a 
factor of 2.4-3 if one adopt [3-2.5. If the flux ratio is calculated 
from observations at 450 and 850 /vm, the experimental ratio k^^q 
I K%s() is equal to 3.6 and 2.6 for the E sample at 10 and 30 K and 
"^450 / '<^850 = 3.6 and 4.9 using the asymptotic expression for the 
dust MAC (Eq|5]l and /? -2 and 2.5 respectively. Hence the flux 
ratio is wrong by a factor of 1-1.4 and 1.4-1.9 if one adopt yS 
equals 2 and 2.5, respectively. The error on the flux ratios is thus 
larger when data at 450 and 1200 //m are used. For the same rea- 
sons it has some consequences on the determination of masses 
from FIR/submm data. At 10 K, the experimental MAC of the 
E sample is /ciioo/j/n ~ 0.45 whereas the value extrapolated using 
the asymptotic expression for k is 0.17, at 850 jjm K^so^m ~ 0.7 
whereas the extrapolated value is 0.45. The error on the mass is 
thus of a factor 3 and 2 at 1200 and 850 yum, respectively. 
The intrinsic MAC for a given material/grain distribution as 
measured in laboratory increases as the temperature increases. 
However, astrophysical observations reveal spatial variations of 
the MAC, and show that the dust emissivity at long wavelength 
seems to increase from the di ffuse medium to dense and cold 
regions by a factor of 3 to 4 dStepnik et alj|2003l : lJuvela et al] 
I2OI U iPlanck Collaboration et al.l 201 la ). Such results seem at 
first sight to contradict the temperature dependence of the MAC 
observed experimentally as a general behaviour in all amor- 
phous samples, but these two experimental and observational 
facts are not necessarily contradictory. To date, the physical 
interpretations associated with these two facts appear different 
and complementary. Experimentally, we measure the variation 
with the temperature of the MAC of a given population of 
specific silicate grains and as a general fact the MAC appears 
correlated with the temperature. Such behaviour is due to 
the intrinsic physical properties of the amorphous state. The 
disorder of the amorphous state induces some additional low 
energy broad band to the MAC that have to be taken into 
account and we cannot consider anymore that the FIR/submm 
MAC in amorphous silicates is only dominated by the long 
wavelength wings of the mid infrared bands. As discuss above 
it can be understood in terms of the TLS model. In contrast, the 
observational spatial variations of the interstellar dust MAC are 
analysed in terms of a spatial evolution of the grain population: 
aggregation, fragmentation, variations of both the size distri- 
bution and the fluffiness or fractal dimension of the grains due 
to the aggregation processes Stepnik et al. (2003). The increase 
of emissivity observed in the core of dense clouds where the 
dust temperature is low probably reflects a variation of the 
grain population, through a change in the composition and the 
morphology of the grains due to the aggregation/fragmentation 
processes and not a variation with the temperature of the 
intrinsic properties of the grain material, as it is observed 
experimentally. Model calculations with the discrete-dipole 
approximation (DDA) can explain this enhancement with the 
coagulation of dust particles into aggregates. In the FIR/mm 
range, such coagulation increases the amplitude of the observed 
emissivity of the dust populatio n without changing its spec- 
tral profile (i.e. spectral index) dOssenkopf & HenningI Il994t 



iKohleretalJ 1201 Ih . However, DDA and similar calculations 
require the knowledge of the intrinsic optical properties of the 
material that constitutes the grains (complex dielectric constant 
or refractive index), and these intrinsic optical properties used 
at the basis of the modelling have to take into account their 
temperature-dependences. 



6.4. Interpretation of /3-T diagrams 

Results from experimen t al work from t his study but also 
from Agi adze et al.l(ll996l) : lMennellaet alJ(ll998h : lBoudet et all 
(2005) have implications on the derivation and interpretation of 
/3-T diagrams that are used to investigate the dust properties in 
various astrophyscial environments. In particular it shows that it 
is difficult, if not impossible, to compare T-/3 diagrams if they 
rely on observations with different spectral coverages. If the 
diagrams are derived from observations whose spectral cover- 
age does not overlap, then their comparison is meaningless be- 
cause it is impossible to conclude anything on the dust prop- 
erties since the same dust population may have very different 
emissivity spectral index in different domains (see Figs. |6]and 
|7]l. However this is almost never the case since data at short 
wavelengths (~ 100 /im) are necessary to constrain the dust 
temperature. If the spectral coverage overlaps, as it is the case 
for example for Herschel and combined Herschel/Planck data 
or Herschel/ground based data, then the derived single emis- 
sivity spectral index, which is averaged over the wavelength, 
is thus different for the two spectral coverages. The more the 
spectral coverages overlap, the more the derived spectral in- 
dexes from the two dataset are comparable. This should be kept 
in mind, in addition to the effect of observational noise on the 
derivation of the dust temperature and emission spectral index 
(Shettv et al. 2009) and to the effec t of temperature variations 
along the line of sight (lJuvela et al.l 1201 1 ). when interpreting 
/3-T diagrams. Nevertheless, there is no doubt from the experi- 
ments presented in this study, but also f r om previous experimen - 
tal studies from Agiad ze et all (Il996h : iMennella et alj (Il998l) : 
Boudet et al. (2005), that an anti-correlation between the silicate 
grain temperature and the spectral index /? intrinsic to the grain 
material does exist (Table [3]l. Hence, even though observational 
bias such as noise and temperature distribution certainly compli- 
cate the data analysis and reduce the accuracy of the derived pa- 
rameters, the /3-T anticorrelation, intrinsic to the grain material, 
exists and certainly contributes to the observed T-yS anticorrela- 
tion. 

Indeed, the experimental MAC variations with temperature are 
in qualitative agreement with astrophysical observations (J3-T 
anti-coiTelation and millimetre excess). However the experi- 
mental /3-T anti-correlation occurs in a wider temperature range 
than the one observed in astrophysics: except for the samples 
with an enstatite stoichiometry, we do not observe significant 
MAC variation within the 10-30 K range. This may have several 
causes. First of all, for all the reasons discussed above, the 
astronomical yS-T anti-correlation can be only partly attributed 
to variation of grains properties with temperature. Secondly, as 
suggested by the experimental results, in agreement with the 
TLS model predictions, the MAC variations are due to disorder 
in the amorphous state (see Sect|5]l. Thus, astronomical silicates 
may present more (or less) pronounced variation of their MAC 
with T and A than the dust analogues studied here, according 
to their micro-structure and composition. The comparison and 
modelling of further astrophysical and experimental observa- 
tions with the same physical model, such as the TLS model. 
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should be useful to discriminate between these points. 



6.5. Perspectives 

From a theoretical point of view, the emission law formulated 
in Eq. [1] for diffuse interstellar dust should be relativized 
since it cannot take into acc ount the dependen ce of the dust 
MAC with the temperature. iMenv et alT (l2007l) proposed an 
alternative physical model for the dust emission, the TLS 
model, based on specific properties of the amorphous state. 
We show that in the disordered structure of the amorphous 
grains, the charge disorder and the density of states of low 
energy levels (which can be seen as a distribution of closed 
rotational-type configurations of atoms or groups of atoms in 
the amorphous structure) lead to absorption processes that can 
dominate, in the FIR/mm wavelength range, over the absorption 
of the mid-infrared resonance wings. This model , which is in 
good agreement with astrophysical observations dParadis et al.l 
l201lh and laboratory-experiments (see Sect. |5] and Meny et 
al. in preparation), is now being improved (Gromov et al. in 
preparation). 

We have presented here a study of Mg-rich amorphous silicate 
grains. However it must be kept in mind that interstellar grains 
are probably more complex than the studied samples. In partic- 
ular they may contain iron in the form of oxide or metallic in- 
clusions or within the silicate network. They may also be coated 
with ices in dense and cold clouds or with carbonaceous matter. 
Grains of carbonaceous matter may also be big enough to radiate 
at thermal equilibrium in the FIR/submm range. These material 
will certainly influence the shape of the dust emission at long 
wavelength. However, the optical properties of ices and carbona- 
ceous matter in the FIR /submm at low temperature are poorly 
known. iMennella et al.l (11998) studied two amorphous carbon 
samples in the range 20-2000 fj.m, from 300 to 24 K. As for the 
silicate grains they found that the MAC decreases when the grain 
temperature decreases and that its spectral shape is slightly mod- 
ified. These experimental results must be taken into account for 
the interpretation of FIR/submm/mm data such as the Herschel 
and Planck data. Carbonaceous dust is not as well constrained as 
the silicate dust and further experimental work have to be done 
to characterise the FIR/submm spectroscopic properties at low 
temperature of this dust component but also of ices, oxides and 
Fe-rich silicate dust. 

7. Conclusion 

In this study we report the characterisation of the far infrared and 
submillimetre optical properties of various interstellar dust ana- 
logues as a function of grain temperature. The studied samples 
are amorphous and crystallised Mg- and Ca-rich nanograins of 
mean composition characteristic of forsterite, enstatite and diop- 
side. For all the amorphous samples, independently of their com- 
position, we observe that: 

1 . The FIR/submm mass absorption efficiency (MAC), or opac- 
ity, decreases when the grain temperature decreases. 

2. The spectral shape of the MAC curve does not follow a sim- 
ple asymptotic behaviour with the wavelength so that it can- 
not be described correctly with a single spectral index. 

3. At a given wavelength, the locally defined spectral index (i.e. 
the slope of the MAC curve) increases when the dust temper- 
ature decreases. 



These variations of the grain MAC, that are not observed for 
crystallised samples, are due to the amorphous nature of the sam- 
ples. However, the detailed spectral shape of the MAC depends 
on the cations/Si ratio of the amorphous materials and is different 
for amorphous materials having olivine and pyroxene stoichiom- 
etry. At short wavelength {A < 500/700jum), the spectral index 
/3 is close to the value of 2 (1.6 < /3 < 2.2) for all grain temper- 
ature and composition whereas at longer wavelengths, /3 < 2 for 
the pyroxene stoichiometric samples (J3 - 0.9 for the Mgo gsSiOs 
sample at 10 K above 800 fj.m) and /3 > 2 for the olivine stoi- 
chiometric samples (J3 - 3.6 for the Mg2.3Si04 sample at 10 K 
above 700 jum). 

These experimental results have important implications for the 
interpretation of astronomical observations. They emphasise the 
fact that adopting a power law characterised by a single spectral 
index independent of the wavelength and the grain temperature 
for the dust MAC is not appropriate and may induce significant 
errors on the derived physical parameters such as the dust tem- 
perature and mass. The study of interstellar dust should be based 
either on the MAC as a function of the wavelength and temper- 
ature rather than on the spectral index. Alternative models such 
as the TLS model that reproduces the variation of the grain opti- 
cal properties with the temperature and wavelength may also be 
considered. Then, new information about the dust, its physical 
and chemical environment and its evolution may be obtained. 
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